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The oscillating scenario of route to Lambda was recently proposed by us [l[ as an alternative 
to a cosmological constant in a explanation of the current accelerating universe. In this scenario 
phantom scalar field conformally coupled to gravity drives the accelerating phase of the universe. 
In our model ACDM appears as a global attractor in the phase space. In this paper we investigate 
observational constraints on this scenario from recent measurements of distant supernovae type 
la, H{z) observational data, CMB R shift and BAO parameter. The Bayesian methods of model 
selection are used in comparison the model with concordance ACDM one as well as with model 
with dynamical dark energy parametrised by linear form. We conclude that ACDM is favoured over 
FRW model with dynamical oscillating dark energy. Our analysis also demonstrate that FRW model 
with oscillating dark energy is favoured over FRW model with decaying dark energy parametrised 
in linear way. 

PACS numbers: 98.80.Es, 98.80.Cq, 95.36.+X 



I. INTRODUCTION 



Observations of distant supernovae type la still consistently suggest that the universe is in a accelerating phase of 
expansion 0, S [3| ■ These confirmations are supported by CMB observations which indicate that universe is almost 
spatially flat [5| and that the amount of matter in the universe calculated from galaxy clustering is not enough to 
account for this flatness [1, 0] ■ These observational facts regarded on the background of standard general relativity 
indicate that about 2/3 of total energy of the universe today being a dark energy with negative pressure which is 
responsible for the current accelerated expansion if the strong energy condition is violated. 

There are many candidates for dark energy description 8, and references therein]. Here we consider dark energy 
in the form of phantom scalar field ij) with the quadratic potential function U{^) for simplicity of presentation. The 
scalar field is conformally coupled to gravity. In our previous work it has been demonstrated that for generic class of 
initial conditions the equation of state parameter WcR = Pcs/PcS approaches —1 value through the damping oscillations 
around this mysterious value. Hence theoretically appeared the possibility to solve the cosmological constant problem 
where the smallness of of cosmological constant does not require fine tuning of model parameters. 

Here we use different astronomical observations to confront the model with the observational data. In this paper 
we use SNIa data and other tests like CMB R shift, BAO and H{z) data obtained from differential ages of galaxies 
0. Bayesian statistics is used to constrain a set of model parameters. In the constraining the model parameters we 
perform combined analysis with CMB R shift parameter as calculated by Wang and Mukherjee for WMAP 3 
The main question addressing in this paper is whether data sets to favour an evolving in oscillatory way dark energy 
model over ACDM one. Using Bayesian framework of model selection we also compare oscillating parametrisation 
with other most popular linear in scale factor a parametrisation. 

Guo, Ohta and Zhang [ll[ developed theoretical method of reconstruction of the quintessence potential directly 
from the effective equation of state parameter w{z) for minimally coupled scalar field. This method can be extended 
to the case of non-minimally coupled scalar field. 
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II. OSCILLATING DARK ENERGY MODEL 

Investigations of different dark energy models [1] are hindered by lack of alternatives to the effective cosmological 
constant model [l^l- The simple step toward more realistic description is that the dark energy might vary in time. 
Usually the form of w{z) is a priori assumption to remove some degeneration problem in analysis of constraints on 
model parameters from observational data. However may happened that assumed form of parametrisation of the dark 
energy equation of state is incompatible with true dynamics which determine w{z) itself. We propose to determine 
corresponding form of w{z) directly from the dynamical behaviour in the vicinity of stable critical point representing 
effective model ACDM. From the dynamical systems methods we know that the system in the phase space can be 
good approximated by its linear part p^ . Then we solve differential equation determining wx{z)- As a result we 
obtain jj] 

wx{z) = -1 + (1 + z)3|ci cos(ln(l + z)) + C2 sin(ln(l + z))| (1) 

for phantom scalar field non- minimally (conformally) coupled to gravity [13] ■ Note that a single scalar field model 
with general Lagrangian L = 9^09^ (/>) will not be able to have w crossing —1 [l^ and to realize that one must 
introduce non-minimal coupling or modification of Einstein gravity. 

We consider conformally coupled phantom scalar field with and given by 



Pip 
Pip 



1 



2H + 3H' 



where dot denotes differentiation with respect to cosmological time. 

From eq.([T]), instead of most popular linear parametrisation, we obtain model with characteristic crossing of wx 
— 1 "phantom divide", thereby the violation of weak energy condition infinite times in the past. 

With the help of formula ([T]) one can simply calculate energy density for dark energy p\ 



PA = PA,o exp{-D2) exp (1 + zy Di sin(ln(l + z)) + D2 cos(ln(l + z)) 



(2) 



where Di = 0.3(Ci + 3C2) and D2 — 0.3(3Ci — C2). It is interesting that some special cases of this dark energy 
parametrisation are explored in probing for dynamics of dark energy [la . flTj . 

Let us consider flat FRW model filled with dark energy with density pA, dust matter (baryonic and dark) and 
radiation. For further analysis of constraints from cosmography it would be useful to write Friedmann first integral 
on H^, where H is the Hubble's parameter 



i? = iJoA/f^A,oexp(-£>2)exp (l + z)3 L>i sin(ln(l + z)) + Da cos(ln(l + z)) + a„,o(l + z)3 + J7^,o(l + ^)^ (3) 



where ^r,o — 0.5 * 10 ^ is fixed and ^Ia,q = 1 ^ ^m.o ~ ^r,Q- Di, D2 and flmfi are free parameters which should be 
fitted from observational data. 



III. CONSTRAINTS FROM SNIA, SDSS, CMB AND H(Z) OBSERVATIONS 

To constrain the unknown values of model parameters we used the set of TV = 192 SNIa data [3, [H, [II- Here 
we based on the standard relation between the apparent magnitude (m) and luminosity distance (di): m — M — 
5 logiQ i^L where M is the absolute magnitude of SNIa, 7W = — 51og^Q i/o + 25 and 13^ = HqcIl- The luminosity 
distance depends on the considered cosmological model and with assumption that fc = is given by = (1 + 

z)c 

Posterior probability for model parameters (after marginalization over nuisance parameter - Hq with the assumption 
that prior probability for this parameter is flat within the interval < 60, 80 >) has the following form 



Pi0\M,D)^ / {7r{e\M) exp [-0.5x1 AS)] )dHo, 



(4) 
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where x|jv(^') = E^i ( ^''"~^''' )^ fJ-f' = m, - M, [if = hXag^^Bu + A^, 71(6*1^/) is the prior probabihty for 
model parameters and d = (r2„j,o, i?2)- Here we assumed fiat prior for model parameters within the interval: 

nm,o e< 0, 1 >, Di e< -1, 1 >, D2 e< -1, 1 >. 

The best fit values for model parameters (the mode of the posterior probability) are the same as the best fit 
values obtained by minimization within the interval for parameters assumed before. Results, i.e. values for model 
parameters obtained via minimization procedure are gathered in Table |T1 
Posterior probabilities for model parameters defined in the following way 



P{nm,o\M,D) = j j P{6\M,D)dDidD2 
PiDi\M,D)= [ [ P{6\M,D)dnm,odD2 



P{D2\M,D)^ j j P{e\M,D)d^,nfidDi 

(5) 

are presented on Figure [TJ The values of the mean for such distributions together with 68% and 95% credible interval 
are gathered in Table HI Two dimensional contour plots representing the 68% and 95% credible interval of the joint 
posterior probability distributions i.e. P{nm.,o,Di\M,D) = J P{e\M,D)dD2, P{n^fi, D2\M, D) = J P{e\M,D)dDi, 
P{Di,D2\M,D) = J P{6\M,D)dVtm.,o are presented on Figure [2 [3] and H respectively. 

We add constraints coming from observational H(z) data (N=9) 0, [13, [Ml. This data based on the differential 
ages (^) of the passively evolving galaxies which allow to estimate the relation H{z) = f = ^T^IH' ^^'^ posterior 
probability for model parameters has the following form 

P{d\M, D)^J {7r{d\M) exp [-OMxInW + XniO))] )dHo, (6) 

where xl,W = Eti (M£Ozp(^)^ 

We also used constraints coming from so called CMB R shift parameter. In this case the posterior probability for 
model parameters has the following form 

Pi9\M, D)^J {tt{0\M) exp [-0.5(x|m(^) + Xh{0) + xlM] )dHo, (7) 

where Xr{0) = (^^^^9^)' and i?*'' = ^O;:^ f^'^^ -^dz, = 1.70 ± 0.03 for zjcc = 1089 

Finally we add constraints coming from the SDSS luminous red galaxies measurement of A parameter (A"'"* = 

0.469 ± 0.017 for za = 0.35) [l^], which is related to the baryon acoustic oscillation peak and defined in the following 

_ 1 A 

way A - V^^.o l^^J^J [~ io 7?(i) ■ 

This parameter was derived with assumption that w{z) is a constant. Due to that using this value to constraints 
varying w{z) lead to systematic errors in the parameter constraints [235. The posterior probability has the following 
form 

P[e\M, D)^ [ {^(e\M) exp [-0.b{xlN{e) + XuiO) + xliO) + x\m)\ )dHo, (8) 



where x^(6') = 



As one can see after inclusion all data to the analysis we obtain the values for Di and D2 parameters which are 
close to zero. Due to that we also consider two models which are special cases of the oscillating dark energy model: 

1. Osc DE 1: C2 = ^ wx{z) = -1 + Ci(l + z)^ cos(ln(l + 0)), 



Hiz) = i?oy f^A,o exp(-0.9Ci) exp |^0.3Ci(l + z)3 ^3cos(ln(l + z)) + sin(ln(l + z)) ] + 17™,o(l + z)^ + 17^,o(l + ^)^ 
where 17^ — 0.5 * 10""* and Ha = 1 ~ ^m.o — o- 
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TABLE I: Values for oscillating DE model parameters obtained via minimization (best fit), values of the mean with the 68% 
and 95% credible intervals obtained from the posterior probability distribution for considered oscillating DE model parameter. 
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FIG. 1: Posterior probability distributions for oscillating DE model parameters. 
2. Osc DE 2: Ci = ^ Wxiz) = -1 + C2(l + z)^ sin(ln(l -|- z)), 



H{z) = Ho^ftAfi exp(0.3C2) exp - 0.3C2(1 + z)^ [ cos(ln(l + ^)) - 3 sin(ln(l + ^))J j + 0„,o(l + z)^ + Or,o(l + z)'^, 
where Clr,o — 0.5 * 10~^ and Qa.o = 1 — ^mfi — ^rfl 
To constrain values of parameters for models defined above we repeat the calculation described before. Results are 



5 




FIG. 2: Contour plots representing the 68% and 95% credible interval of the joint posterior probability distribution for Qm.,o 
and Di oscillating DE model parameters. 



gathered in Table |TT] and IIIII respectively. Posterior probabilities are presented on Figure [5] and [7] respectively. Two 
dimensional contour plots in the {flm,o, Ci) plane are presented on Figure [6] and [8] respectively. 
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TABLE II: Values for oscillating DE 1 model parameters obtained via minimization (best fit), values of the mean with the 
68% and 95% credible intervals obtained from the posterior probability distribution for considered Osc DE 1 model parameter. 



The wxiz), and — functions together with 68% credible interval for considered models ( calculated for the 

^ PA,0 Pm ° ^ 

mean of the posterior distributions for the model parameters which are gathered in TableHlUl l IIIII in the SN+H+R+A 
case) are presented on Figure El [TOl HH respectively. 
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FIG. 3: Contour plots representing the 68% and 95% credible interval of the joint posterior probability distribution for Qm.,o 
and D2 oscillating DE model parameters. 
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TABLE III: Values for oscillating DE 2 model parameters obtained via mirnmization (best fit), values of the mean with the 
68% and 95% credible intervals obtained from the posterior probability distribution for considered Osc DE 2 model parameter. 



Finally we made a comparison of Oscillating DE Models, ACDM model and model with linear in scale factor 
parametrisation of w. w{a) ^ wq + wi{l — a). Analysis was made in the Bayesian framework. Here the best model is 
this one which has the largest value of the posterior probability. It is convenient to use the posterior odds in analysis, 
which in the case when no model is favoured a priori is reduced to so called Bayes Factor Bij (the ratio of the evidence 
for models indexed by i and j) [13, [IE]- This quantity can be interpreted as a strength of evidence against worse model 
with respect to the better one: < \2\n Bij\ < 2-not worth more than a bare mention, 2 < 12 In | < 6 ~ positive, 
6 < 1 2 In By I < 10 - strong, and \2\n Bij\ > 10 - very strong. Here we used BIC quantity [2^] as an approximation 
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FIG. 4: Contour plots representing the 68% and 95% credible interval of the joint posterior probability distribution for Di and 
D2 oscillating DE model parameters. 

to the minus twice logarithm of the evidence, which is defined in the following way: 

BIC = -2\nC + d\nN, 

where C is the maximum of the likelihood function, d is the number of model parameter and N is the number of data. 
Values of Bayes Factor (calculated with respect to ACDM model) are gathered in Table IIVI 



Model 


21nB 


ACDM 





Osc DE 


8.75 


Osc DE f 


3.37 


Osc DE 2 


3.37 


Linear parametrisation 


5.76 



TABLE IV: Twice logarithm of the Bayes Factor 

As one can conclude ACDM model is the best one from the set of models considered in this paper. Evidence in 
favour this model is strong when comparing with the Osc DE model and positive in the other cases. There is positive 
evidence in favour model with Linear parametrisation over the Osc DE model. Bayes Factor computed for Osc DE 1 
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FIG. 5: Posterior probability distributions for Osc DE 1 model parameters. 



and Osc DE 2 models is close to 1 which indicate that the information coming from the data (used in analysis) are 
not enough to favour one of this model over another. In this situation calculation the Bayesian evidence by numerical 
integration could give better results. Finally Osc DE 1 and Osc DE 2 are favoured over the Osc DE Model and over 
model with linear in a parametrisation of w. 



IV. CONCLUSIONS 



In this paper we have placed constraints on a parametrised dark energy model [l| using the SNIa data sets, 
observational II(z) data, the size of the baryonic acoustic oscillation peak from SDSS and the shift parameter from 
the CMB observations. We study possibility that phantom dark energy is oscillating rather than decaying to A. Such 
a scenario opens the possibility of the non-minimal coupling to gravity for phantom scalar field. Combining four data 
bases (SNIa, H(z), CMB, SDSS) we obtain constraints on the oscillating dark energy model parameters (Om, Di, D2) 
and compare this model with ACDM model and with model with linear in a parametrisation of w in the Bayesian 
framework. It is found that special cases of oscillating phantom dark energy model ( called Osc DE 1 and Osc DE 
2 in this paper) are favoured over the model with linear in a parametrisation of w. Cosmological constant case still 
remains as the best one from the set of considered models. 
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FIG. 6: Contour plots representing the 68% and 95% credible interval of the joint posterior probability distribution for Qm,o 
and Ci Osc DE 1 model parameters. 
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FIG. 8: Contour plots representing the 68% and 95% credible interval of the joint posterior probability distribution for Qm.,o 
and C2 Osc DE 2 model parameters. 




FIG. 9: wx{z) for the Osc DE, Osc DE 1, Osc DE 2 Model (black line) together with 68% credible intervals (dashed lines) 
(calculated for the mean and 68% credible interval of the posterior distributions for the model parameters) and for the ACDM 
model (red line). 
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FIG. 10: for the Osc DE, Osc DE 1, Osc DE 2 Model (black line) together with 68% credible intervals (dashed Unes) 

(calculated for the mean and 68% credible interval of the posterior distributions for the model parameters). 




FIG. 11: ^ for the Osc DE, Osc DE 1, Osc DE 2 Model (black line) together with 68% credible intervals (dashed lines) 
(calculated for the mean and 68% credible interval of the posterior distributions for the model parameters). 



